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Abstract 

Background: Ionizing radiation is genotoxic to cells. Healtliy tissue toxicity in patients and radiation resistance in tumors 
present common clinical challenges in delivering effective radiation therapies. Radiation response is a complex, polygenic 
trait with unknown genetic determinants. The Drosophila Genetic Reference Panel (DGRP) provides a model to investigate 
the genetics of natural variation for sensitivity to radiation. 

Methods and Findings: Rad'\at\or\ response was quantified in 154 inbred DGRP lines, among which 92 radiosensitive lines 
and 62 radioresistant lines were classified as controls and cases, respectively. A case-control genome-wide association 
screen for radioresistance was performed. There are 32 single nucleotide polymorphisms (SNPs) associated with radio 
resistance at a nominal p<10~^; all had modest effect sizes and were common variants with the minor allele frequency > 
5%. All the genes implicated by those SNP hits were novel, many without a known role in radiation resistance and some 
with unknown function. Variants in known DNA damage and repair genes associated with radiation response were below 
the significance threshold of p<10~^ and were not present among the significant hits. No SNP met the genome-wide 
significance threshold (p = 1.49x10"^), indicating a necessity for a larger sample size. 

Conclusions: Several genes not previously associated with variation in radiation resistance were identified. These genes, 
especially the ones with human homologs, form the basis for exploring new pathways involved in radiation resistance in 
novel functional studies. An improved DGRP model with a sample size of at least 265 lines and ideally up to 793 lines is 
recommended for future studies of complex traits. 
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Introduction 

Ionizing radiation is genotoxic to cells. The acute effects of 
whole body irradiation depend on several factors, among which 
cell type and total dose are most prominent. Whole body exposure 
to radiation is implicated in gastointestinal syndrome and 
hematopoietic syndrome, and for doses higher than 30 Gy, 
cardiovascular collapse, central nervous system damage and death 
within 24—72 hours [1]. Direct cellular response to radiation 
involves formation of complex DNA double strand breaks, leading 
to mutations and genomic instability that may cause cancer [2]. 
The progeny of irradiated human and murine cells have been 
demonstrated to develop chromosomal abnormalities or other 
mutations after multiple generations, and these genomic changes 
are the same as those observed in human tumors [3] . Thus, there 
are strong genetic components involved in radiation response. 
Interestingly, there is no evidence for radiation-induced germline 



mutations or heritable genetic diseases in children of irradiated 
parents [4,5,6]. In comparison, somatic cell radiosensitivity has 
been established to be a heritable trait in humans [7,8]. 

It has been estimated that approximately 80% of inter-patient 
variation in normal tissue toxicity following radiation exposure is 
likely to be genetic, whereas only 20% of the variation results from 
stochastic events associated with the random nature of radiation- 
induced cell killing in addition to random variations in dosimetry 
and dose delivery [9] . Studies addressing this have mainly focused 
on candidate genes, however they have failed to demonstrate 
vmequivocal links between genotype and radiation toxicity [10]. 

Radiotoxicity, also known as radiosensitivity, is thought to be a 
polygenic trait involving interactions of multiple loci involved in 
different cellular pathways, most of which confer a relatively small 
risk [10,11]. To date, three human genome-wide association 
(GWA) studies have been conducted on radiation response. 
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leading to the discovery of two new genes, FSHR and PRDMl, 
involved in radiotoxicity [12]. However, human GWA studies 
suffer from a number of limitations, the chief one being the 

inability to identify causal variants [13,14]. 

The Drosophila Genetic Resource Panel (DORP) is a powerful 
community resource for interrogating heritable, natural variation 
and linking complex traits to underlying genotypes [15]. Spon- 
sored by the National Human Genome Research Institute, this 
open platform contains 192 fuUy sequenced inbred Drosophila 
strains harboring nearly 5 million single nucleotide polymorphisms 
(SNPs). As a discover^' platform, the DGRP offers compelling 
advantages for understanding how genes specify complex traits. 
First, all sequence data, SNP calls, stocks and associated web- 
based tools are publicly available to the research community. 
Second, unlimited sampling of these highly inbred strains can 
dramatically magnify statistical power. As experimental surrogates 
for individual variation, DGRP stocks collectively deliver a much 
higher statistical power compared to outbred individuals [16]. 
Third, nucleotide level resolution delivers sequence variants that 
are more likely to be causative [15], which distinguishes it from 
human GWA studies that deUvers SNPs that are usually linked to 
an unknown culprit. A final advantage afforded by the DGRP 
draws on the power of the model itself As a genetic system, 
Drosophila offers opportunities for follow up examination of 
candidate variants. Taken together, these capabilities facilitate 
unbiased discovery and functional examination of natural 
sequence variants that determine traits of interest. 

Here, we conducted GWA mapping of radiation resistance in 
the DGRP lines by assaying survival of adult males at a fixed 
radiation dose. 

Materials and Methods 

Drosophila stocks 

We used 154 inbred lines of the DGRP [15], available from the 
Bloomington Drosophila Stock Center. Flies were reared at a 
controlled density on cornmeal-molasses-agar medium at 25°C, 
60-75% relative humidity and a 12-h light-dark cycle. 

Experimental setup 

The radiation source was a commercial Cs' " irradiator (J. L. 
Shepherd, Model Mark I-68A, Serial Number 1158, San 
Fernando, CA) emitting 662 keV gamma photons. The irradiator 
was calibrated using a PTW N31010 ionization chamber (PTW - 
New York Corporation, Hicksville, NY). The air-kerma calibra- 
tion coefficient for Cs ' photons was obtained from the 
Accredited Dosimetry Calibration Laboratory (ADCL) at the 
University of Wisconsin. In addition to the absolute dose rate 
calibration at a reference point, a relative isodose distribution 
within the irradiator was obtained using radiochromic films 
(Gafchromic EBT3, International Speciality Products, Wayne, 
P\J). The flies were irradiated in cylindrical vials with an inner 
diameter of 2.3 cm, where the lower part of the vial typically had 
2.5 cm of food and the flies were able to fly within 5 cm above 
that. The vials were placed at the bottom of a rotating plate and 
the plate was kept rotating the whole time during irradiation. The 
average dose rate in the vial was 4.85 Gy/min, resulting in a total 
dose of 1 382 Gy over the 4 hours and 45 minutes of continuous 
exposure. 

Radiation assay 

We developed and optimized a radiation toxicity assay using the 
Canton S strain, controlling for the effects of age (7-11 days old), 
sex (males only), number of flies per trial (n = 50) and survival 



phenotype endpoint (24 hours post-irradiation). We also found no 
noticeable difference in vial temperature before and after 
irradiation, eliminating resistance for heat stress as a potential 

confounder. 

We used the assay in a pilot experiment to test for its ability to 
induce variation in response in 10 randomly chosen DGRP lines. 
We then measured the radiation response of the remainder 144 
lines, and collected data for two replicates per line. For variable 
Hnes, we performed up to 7 replicates per hne. Survival was 
defined as the abihty of males to fly 24 hours post-irradiation, and 
survival time was expressed in percentage as the mean of two 
trials. 

Phenotypic stability 

For temporal phenotypic stability experiments, 12 highly 
resistant DGRP lines were maintained for five to fourteen months 
under normal rearing conditions, and were reassayed for radiation 
resistance using the identical assay conditions. 

Extreme phenotype line crosses 

Reciprocal matings were set up between completely sensitive 
line RAL-28 and extremely resistant line RAL-69, and Fls were 
selfed to generate F2s. Both Fls and F2s were aged 7-1 1 days and 
their racUation response measured under identical assay condi- 
tions. 

Genome-wide association analyses 

We treated radiation response phenotype as a binary outcome - 
resistant and sensitive - since 60% of the data had numerical 
values of zero (Figure 1). A subset of resistant lines was highly 
variable (see Results). However, variability did not affect our data 
analysis because variable lines exhibited a mean resistant 
phenotype after several replications. Similarly, although the 
temporal phenotypic stability analysis also showed variability in 
highly resistant lines, no highly resistant line became sensitive. 
Hence, variabihty in the resistance phenotype did not alter the 
number of cases and controls for association analysis. 

There were a total of 5,066,519 SNPs in the DGRP freeze 1 
dataset. Genotype data were cleaned by the following criteria: 
genotype missingness <15%, heterozygous haploids and genotype 
call rate >10% and minor allele frequency >1%. There were 
2,035,449 SNPs filtered out by quality controls. 

Potential confounding factors - Wolhachia infection status, 
population stratification, and cryptic relatedness — were consid- 
ered. Wolhachia infection status had no significant effect (Table 
S3). To examine whether population structure is an influential 
confounder, we first derived the top five principal components 
using all the 2,035,449 SNPs by GCTA [17]; then we tested 
association between radio resistance and each principal compo- 
nent. The results suggested population structure not a significant 
confounder (Table SI). To thoroughly control for population 
stratification and cryptic relatedness, we employed a linear mixed 
model that uses the whole-genome data to estimate the genetic 
relationship matrix, as well as the top five principal components as 
covariates [18]. In sum, the association of radiation response in 
154 lines vidth 3,030,570 SNPs was examined by the likelihood 
ratio test [19] fitting a linear mixed model using GEMMA [20]. 

To determine the threshold for genome-wide significance, we 
calculated the number of haplotype blocks in the DGRP genome. 
The haplotype block is defined as a window of SNPs with the 
outer-most marker required to be in strong linkage disequilibrium 
(LD) with an upper limit of 90% confidence interval exceeding 
0.98 and a lower limit of 90% confidence interval exceeding 0.7 
[21]. The calculation was performed using PLINK [22]. Note that 
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Figure 1. Variation in radiophenotype among 154 DGRP lines. 50 males from each line were aged for 7-1 1 days, irradiated with 1382 Gy and 
the number of survivors scored 24 hours post-irradiation. The data shown is the mean of two independent trials; the error bars represent standard 
deviation. The 62 resistant and 92 sensitive lines are indicated by blue lines. 
doi:1 0.1 371 /journal.pone.01 04858.g001 



the pairwise LD was only calculated for SNPs within 500 kb; thus 
the number of LD blocks obtained would be an upper limit. The 
narrow sense heritability from additive genetic effects was 
estimated by a hability threshold model using the whole-genome 
variation [23]. We also estimated heritability by fitting traditional 
linear mixed models modelling line effects. 

Results 

Inbred strains exhibit natural variation in radiation 
response 

To characterize natural variation in radiation response in 154 
DGRP lines, we selected a fixed dose of 1382 Gy and Canton S 
strain as a test case. In preliminary experiments we found that 
male response was less variable than female response, hence we 
used males only for our study. 

To examine the suitability of our radioassay for the DGRP lines, 
we selected a panel of 10 DGRP lines at random as well as 3 
common Drosophila lab strains with distinct genetic backgrounds 
(Canton S, yw and w'"^). The variation in the DGRP lines 
covered almost the entire spectrum of radiotoxicity (0-95%) with 
very low standard deviation, while the observed variation in the 
lab strains ranged from 2-28% (Fig. SI in File SI). 

We selected two DGRP lines with extreme phenotypes from our 
pilot panel, RAL-69 and RAL-28, set up reciprocal crosses and 
scored the radiotoxicity of Fl and F2 males. The parental strains 
exhibited a similar radiation response as before (0% survival for 
RAL-28, 95% survival for RAL-69; Fig. 2). The FIs from both 
crosses showed near-sensitivity of radiation response, while the F2s 
from both crosses exhibited at least three-fold higher survival than 
the FIs (^ = 8.06x10"'' by ANOVA). Fitting the parental strains 
and FIs' survival to different genetic models identified the most 
parsimonious model to be a recessive model for radioresistance 
(Table S4.1). The proportions of survival in F2s were 0.36, 0.14, 
0.20, and 0.24, which is close to a recessive proportion of 0.25. 



Therefore, the results suggest a recessive mode of inheritance for 
radioresistance for line RAL-69. 

Because our assay scored survivors 24 h post-irradiation and 
since animals that were alive at the 24 h endpoint showed normal 
life spans, we asked whether the assay was simply finding some 
strains that took longer to die. To address this, we irradiated two 
highly resistant strains, RAL-91 and RAL-142, and scored the 
survivors for up to 10 d (Fig. S2 in File SI). After 10 d, RAL-142 
showed 93% survival, and RAL-91 70%. Also, the flies were 
normal in their behavior and no visible mutations were observed. 

We then screened the remaining DGRP lines for survival 24 h 
post-irradiation, and found extensive phenotypic and genetic 
variation in radiation response (Figure 1, Table S2). We found 
that the lines fell into two distinct groups: one group of 92 lines 
with no survivors at all, another group of 62 lines with survivors. 
We designated the former group as sensitive and the latter as 
resistant. Among the resistant group, there was extensive 
phenotypic variation (ranging from 1% to 98%, Table S2) and a 
subset of lines were highly variable (coefficient of variation >25%). 
To test the stability of the radiation resistance phenotype, we 
selected 1 2 highly resistant lines and retested them after a period of 
five to fourteen months (Table S5). The Spearman's rank 
correlation coefiicient between the two measures was 0.56 with 
a marginally significant p-value of 0.06. The measures in the 
second test were systemically lower than those in the first test for 
unclear reasons (/^ = 4.78 x 10 by a paired t-test); however, all the 
lines maintained their radioresistance. For the resistant lines that 
were highly variable (RAL-149, -237 and -378), we performed up 
to 7 replications for 1 2 lines but they still remained variable (data 
not shown). We speculated that these lines were inherendy 
variable. 

The DGRP lines vary in Wolhachia infection status [15], but 
this parameter had no significant effect on the overall DGRP 
radiation response (Tables S3 and S6). 
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Figure 2. DGRP radloresistance is heritable. Reciprocal crosses between a completely sensitive RAL-28 and a highly resistant RAL-69 lines were 
set up to generate F1, which were then selfed to produce F2. 50 males from F1 and F2 of each cross were scored for survival after 1382 Gy irradiation. 
The data shown represents the mean of two independent trials. 
doi:1 0.1 371 /journal.pone.01 04858.g002 



SNPs associated with radiation resistance 

To identify genes that harbor alleles conferring inter-individual 
differences in radiation resistance, we performed a case-control 
GWA analysis of the phenotype using SNPs from the DGRP 
freeze 1 sequencing data [15]. Since our extreme phenotype line 
crosses demonstrated that radloresistance is a recessive trait 
(Fig. 2), we defined resistant lines as cases and sensitive ones as 
controls. There were 3,030,570 SNPs that met the quality control 
thresholds were tested for association with radloresistance using a 
linear mixed model. 

The quantHe-quantUe plot of p values showed no systematic 
deviation from the null distribution (Fig. 3A). A total of 334,729 
LD blocks were defined, which led to a genome-wide significance 
threshold of 1.49x10"'. The lowest p value in our screen was 
3.66x10"', with none of the SNPs exceeding the genome-wide 
significance threshold. We estimated the sample size needed to 
exceed ^ = 1 .49 x 10"' with a power of 0.8 for the top 32 SNPs, 
and found the required sample size to range from 265 (for 
3R_10180595) to 793 (for 2R_13959724). As the DGRP consists 
of lines derived from a natural population, the prevalence of 
radloresistance was set to be 0.4 (62/154). The heritabUity using 
the whole-genome variation was estimated to be 0.067; in contrast, 
by linear mixed models accounting for line effects, the heritabUity 
was estimated to be >0.8 (Table S4). 

A total of 32 SNPs were associated with radiation resistance at 
p<lO ^, among which two were at ^<10 ^ (Table S7; Fig. 3B). 
The majority of the 32 SNPs were common variants, having minor 
allele frequency of at least 0.14. Their effect sizes (odds ratios) were 
similar, ranging from 1.19 to 1.44 (Fig. 4). 

According to site class, of the 32 SNPs, 15 were intronic, 8 fell in 
coding regions (two non-synonymous, 6 synonymous), two 
occurred in 3' UTR, 6 were intergenic and one was located in a 
gene desert. Excluding the one in the gene desert, the remaining 
SNPs implicated 24 genes, 9 of which have human homologs. AU 
24 genes are novel candidates in radiation resistance (Table S8). 



Pathway analysis using the DAVID software (http://david.abcc. 
ncifcrf gov/) showed that the 32 SNPs were not enriched in any 
known biological process or pathways involved in radiation 
resistance. Furthermore, the 9 human homologs also were not 
enriched in any known biological pathways. 

To determine whether canonical DNA damage and repair gene 
variants were present or enriched among the significant hits in our 
screen, we considered a comprehensive list of 102 genes involved 
in DNA damage response, of which 97 had human homologs, 
based on the compilation done by Wood et al. and Lange et al. 
[24,25,26]. The DGRP set harbored a total of 10,916 variants 
(excluding intergenic SNPs) in these genes (Table S9). We found 
that none of the DNA damage and repair gene variants were 
present among the significant SNPs (Fig. 5). Table SIO lists the 
lowest p value variant of each of the 102 DNA damage and repair 
genes. 

Discussion 

Genes involved in radiation toxicity/ resistance have largely 
been identified using either cultured cells or laboratory-engineered 
single gene mutation studies in animal models. GWA studies have 
emerged as a powerful approach to circumvent the limitations of 
cultured cells and single gene mutation studies, uncovering novel 
loci involved in a number of complex traits that play a role in 
whole organism [20]. We therefore leveraged the publicly 
available DGRP resource harboring natural variation in Drosoph- 
ila to discover novel variants involved in radiation resistance by 
using a whole animal radioassay. Although ionizing radiation is 
not expected to be strong in contemporaneous environments, and 
is not an important agent of natural selection in the wHd, it is an 
important consideration for manned space expeditions. 

This is the first initiative to approach determinants of radiation 
resistance in whole animals. As a heritable and complex trait, post- 
irradiation survival is clearly an intrigiaing phenotype. We 
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Figure 3. Association analyses of radiation resistance among 154 DGRP lines. (A) Quantile-quantile plot. Tine red line indicates the 
expected and the black line the observed p values. (B) Manhattan plot of p values. The red dashed line indicates p<10^^. 
doi:1 0.1 371 /journal.pone.01 04858.g003 



considered several experimental phenotypes with aim to devise an 
experimental design in which feasibility was a pivotal consider- 
ation. For example, we counted the number of offsprings of 
radiation resistant and sensitive flies as an endpoint, and found no 
significant difference (data not shown). Next, we considered 
estabhshing the LD50 values for all 192 strains but quickly found 
that this would not be practical. We also considered larval 
irraidation, where the endpoint would be larval survival or death; 
however, it is well documented that different strains can differ with 
respect to developmental pace. Furthermore, radiation sensitivity 
is highly variable at different stages in development. Therefore, by 
challenging mature adults we were able to focus the study on our 
trait of interest (radio-toxicity) without introducing confounding 
variation derived from differences in developmental rates. 

Our screen is the largest effort to date utilizing the natural 
variation in a model organism to examine the genetic basis of 
radiation resistance as a complex trait. Our phenotype was 
heritable and robust, measuring true survival rather than 
variations in time-to-death (Fig. 2 and S2). For practical reasons, 
our protocol included just two replicates per line. Therefore, 
heritabUity estimation using traditional linear mbced from this data 
are severely limited. We speculate the line effects and random 
errors are not identifiable from each other given two replicates per 



line. The minimal heritabHity estimate from the whole-genome 
variation was consistent with the GWA screen results of lacking 
significant signals; it might suggest a bulk of non-additive genetic 
effects on radioresistance [27]. 

We discovered 24 novel Drosophila genes (of which 9 have 
human homologs), encompassing a total of 3 1 SNPs, likely to be 
involved in radiation resistance, as well as a gene desert SNP that 
might play a regulatory role. None of the 10,916 variants in 102 
canonical DNA damage and repair genes were present among the 
32 significant hits. This finding is consistent with the failure of 49 
candidate-gene studies covering 3,144 SNPs in 1,494 human genes 
to discover variants associated with radiation sensitivity [10], while 
only 3 human GWA studies on radiotoxicity have identified 
variants in two novel genes, FSHR and PRDMl [12]. Therefore, 
our results support the hypothesis that genetic determinants of 
whole animal radio-toxicity could radically differ from determi- 
nants accessed through cultured cells or single gene mutant animal 
studies. 

Among the genes represented by the significant SNP hits, four 
contain multiple SNPs (NK7.1, lack, pnr, CG14621). Future 
studies will aim to ascertain whether, biologically, they have been 
truly selected or are simply in linkage disequilibrium with causative 
variation in other loci. The genes we identified are not within the 
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Figure 4. Odds ratio histogram of top 32 SNPs at p<10~^. 

doi:1 0.1 371/journal.pone.01 04858.g004 

QTL regions Q,l, 0,2, Q3, MTl and MT2 identified by Gomez et 
al [28]. 

SNPs identified in human GWA studies tend to have a modest 
effect size ranging from 1 to 3, with the majority having odds ratios 
in the range of 1 to 1.1 [29]. The effect sizes of significant hits in 



our screen were also modest, ranging from 1.19 to 1.44. Our 
results demonstrate that even an inbred, model system such as the 
DORP, in which each loci is homozygosed and no SNP is left 
unsequenced, does not provide us with the required sensitivity and 
specificity for individualized prognosis. This suggests that the goal 
of personalized medicine is very, very distant — if achievable at all 
- and that precision medicine should be the aim of future GWA 
studies [29,30]. 

Our screen did not yield any SNP that passed the genome-wide 
significance threshold ot p value< 1.49x10 ' based on the LD 
structure of the Drosophila genome, suggesting that the DGRP 
resource is underpowered for this trait. We would need a 
minimum sample size of 265 for at least one among the 32 
significant hits and a maximum of 793 lines for all 32 SNPs to pass 
the genome-wide significance threshold. However, the DGRP 
resource is adequately powered for complex traits such as 
oxidative stress resistance and tunicamycin-induced ER stress 
resistance [31,32]. This suggests that certain complex traits, such 
as radiation resistance, require higher DGRP sample size than is 
available. An important objective for future complex trait research 
is to anticipate features of dichotomous traits under study that wiU 
allow us to predict power requirements early on in the study, prior 
to a full-scale phenotypic analysis. 

It is highly unlikely that ionizing radiation was a source of 
selective pressure during fly evolution and, therefore, we do not 
propose that our work is accessing a 'defense against radiation'. 
We stress that the same holds true for patients exposed to this same 
stressor in the clinic, despite diflFerential sensitivities that are, in 
part, genetic. Like most gene-directed differences in perturbation 
responses, ionizing radiation is probably a surrogate for other 
environmental challenges that mold complex traits. In fact, an 
important motivation for this study was to extract informative 
clues that could advance our understanding of the underlying 
radiation biology. We have not investigated radiation-stable 
proteins in the present study and certainly it is possible that 




Figure 5. Variants in DNA damage and repair genes are not among the top associated SNPs. A scatterplot of p values for all 1 0,91 6 SNPs 
representing a comprehensive set of 102 DNA damage and repair genes (dots), along with the p values of top 32 SNPs (crosses). 
doi:1 0.1 371/journal.pone.01 04858.g005 
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protein stability plays a part in radiation response. Whatever the 
ultimate mechanism turns out to be, the underlying architecture 
governing these differences are distinct when accessed in whole 
animals versus cultured cells. 
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Table SI Association between radioresistance and the top five 
principal components derived from the whole genome SNP data 
of 154 DGRP lines. 

(DOCX) 

Table S2 Raw and mean survival values of 154 DGRP lines. 
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Table S3 Wolhachia infection status and mean radiation 

response values of 154 DGRP lines. 

(DOCX) 

Table S4 HeritabUity estimation of radioresistance. Table 
S4.1. Inferring the genetic model of radioresistance by reciprocal 
cross of the sensitive line RAL-28 and the resistant line RAL-69. 

(DOCX) 

Table S5 Temporal phenotypic stability study of 12 highly 
resistant DGRP lines. 

(DOCX) 

Table S6 Wolhachia infection status has no signficant effect on 
the overall DGRP radiation response. 

(DOCX) 

Table S7 GWA analysis results. 
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Table S8 Functional annotation of top candidate genes. 
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